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Summary 

A full-engine simulation of the three- 
dimensional flow in the GE90-94B high bypass 
ratio turbofan engine has been achieved. It would 
take less than 1 1 hr of wall clock time if starting 
from scratch through the exploitation of parallel 
processing. The simulation of the compressor 
components, the cooled high-pressure turbine, and 
the low-pressure turbine was performed using the 
APNASA turbomachinery flow code. The com- 
bustor flow and chemistry were simulated using 
the National Combustor Code (NCC). The engine 
simulation matches the engine thermodynamic 
cycle for a sea-level takeoff condition. The simula- 
tion is started at the inlet of the fan and progresses 
downstream. Comparisons with the cycle point are 
presented. A detailed look at the blockage in the 
turbomachinery is presented as one measure to 
assess and view the solution and the multistage 
interaction effects. 


Introduction 

One of the goals of the Numerical Propulsion 
System Simulation (NPSS) Program at the NASA 
Glenn Research Center has been to demonstrate a 
high-fidelity three-dimensional turbofan engine 
simulation. This simulation will support the mult- 
dimensional, multifidelity, multidiscipline concept 
of the design and analysis of propulsion 
systems for the future. This paper describes the 


current status of one major part of that goal: the 
complete turbofan engine simulation using an 
advanced three-dimensional Navier-Stokes turbo- 
machinery solver, APNASA, coupled with the 
National Combustion Code (NCC). 

A production engine has been chosen for this 
demonstration: the GE90 turbofan engine shown 
in figure 1. A sea-level, Mach 0.25, takeoff condi- 
tion has been used for the simulation. The main 
reason is that detailed cooling flows for the turbine 
are well known at takeoff since this represents the 
cooled turbine design condition. Since the cooling 
flow represents a significant amount of the re- 
quired boundary condition information, it was felt 
this was a good point for the simulation. It also 
represents a condition where there are the highest 
temperatures and most stress in the engine, and 
is therefore a practical point to gain further 
understanding. 

The GE90 development program included 
component testing of all the turbomachinery as 
well as the combustor. The full engine simulation 
effort has taken advantage of this. All the turbo- 
machinery components have been analyzed and 
compared to component test data to validate and 
calibrate the approach. These efforts have been 
presented by Adamczyk (ref. 2), Turner et al. 
(ref. 11), and Turner (ref. 12). Also presented by 
Turner et al. (ref. 11) was the coupled high- and 
low-pressure turbine run at engine conditions. 
Turner (ref. 13) has presented an uncoupled fan 
and booster with a coupled high-pressure com- 
pressor (HPC) combustor, high-pressure turbine 
(HPT), and low-pressure turbine (LPT). This 
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Figure 1 . — GE90 turbofan engine. 


paper expands on that work by coupling all the 
components. 

One other important aspect of this effort is to 
reduce the wall clock time to obtain a simulation 
by exploiting parallel processing as much as pos- 
sible. Reduced solution times are necessary to 
efficiently explore the physics of a problem, to 
validate a code, and most importantly, to impact 
design. 

The approach used is presented in the next 
section of this paper that describes the GE90 en- 
gine, the turbomachinery solver, the combustor 
solver, and the coupling strategy. The aerody- 
namic blockage is derived which is a catchall for 
three-dimensional and multistage effects. It is a 
useful quantity for assessing and visualizing the 
solution. The results section briefly discusses the 
simulation compared with the cycle, the solution 
timings, and the aerodynamic blockage effects. 

Symbols 

A annular area 
h enthalpy 
P pressure 
r radius 
R gas constant 
T temperature 
V velocity 
m mass flow rate 
p density 

co shaft rotational speed 
y specific heat ratio 

Subscripts and superscripts: 

h hub 
t tip 

T total or stagnation conditions 
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0 tangential 
z axial 

mass-averaged 

mass-averaged 
derived quantity 

Approach 

The overall approach is to simulate the turbo- 
machinery components using the APNASA code 
that incorporates the average passage model. The 
compression system is currently broken into three 
domains: the fan, booster, and HPC. Profiles are 
passed from one of these domains to the next. 
Both the HPT and LPT are simulated in one do- 
main. The combustion simulation uses the parallel, 
unstructured code (NCC). Coupling between these 
codes is across common interface planes. 


The GE90 Turbofan Engine 

The GE90-94B is a production engine offered 
on the Boeing 777-200ER as shown in a cutaway 
schematic in figure 1 . It is a 94 000 lb thrust ver- 
sion of the GE90 with a bypass ratio of 8.4. The 
fan is 120 in. in diameter comprising 22 composite 
wide-chord blades. The fan outlet guide vane 
(OGV) has several types with different camber to 
guide the bypass flow around the pylon. Only the 
nominal type is modeled in this simulation. 

The booster, or low-pressure compressor, 
consists of three stages and rotates on a common 
shaft with the fan and LPT. The three stages com- 
prise seven blade rows. A fan frame strut separates 
the booster from the HPC. The 10-stage HPC has a 
23 to 1 pressure ratio and has been redesigned 
using three-dimensional aerodynamic (3-D Aero) 
technology. The original GE90 HPC was a scale 
of the HPC developed for the NASA/GE Energy 
Efficient Engine (E 3 ) program. The HPC is on a 
common shaft with the HPT. 

The combustor is a dual dome annular design 
for reduced NO x emission levels and reduced un- 
bumed hydrocarbon, carbon monoxide, and smoke 
levels. There are 30 pairs of fuel nozzles around 
the annulus. 

The HPT has two stages comprising four blade 
rows. There is a turbine midframe strut separating 
the HPT and LPT. It diffuses the flow through a 


high-angle outer diameter casing to a high diame- 
ter for an improved efficiency LPT. The LPT has 
6 stages comprising 12 blade rows. A turbine rear 
frame strut follows this. 

The simulation consists of 49 blade rows of 
turbomachinery and a 24° sector of the combustor. 
The 49 blade rows include the fan, bypass splitter, 
nominal OGV (the pylon and different OGV types 
are not modeled), the 3 -stage booster (7 blade 
rows), the fan frame strut, the 10-stage HPC 
(21 blade rows), the 2-stage HPT (4 blade rows), 
the turbine midframe strut, the 6-stage LPT 
(12 blade rows), and the turbine rear frame strut. 
Although there are 30 pairs of fuel nozzles, actual 
periodicity of the geometry requires the modeling 
of two pairs of fuel nozzles, or a 24° sector. 
Unique cooling arrangements make one pair of 
fuel nozzles a-periodic. 

Turbomachinery Simulation 

The average passage approach is incorporated 
in the code APNASA as discussed by 
Adamczyk et al. (ref. 1). Adamczyk (ref. 2) in 
another paper also discusses the use of this tool to 
support turbomachinery aerodynamic design. 
Kirtley et al. (ref. 5) generalized the closure and 
allowed grids with less shear to be used. Some of 
the numerical details used in APNASA are 

1. Four-stage Runge-Kutta explicit three- 
dimensional RANS solver 

2. Local time steps 

3. Implicit residual smoothing 

4. Implicit k-E turbulence model 

5. Two levels of parallel using MPI message 
passing 

Some of the physical models APNASA can handle 
are 

1. Multistage effects modeled using the aver- 
age passage formulation of Adamczyk and 
generalized for nonpure H-grids 

2. Wall functions 

3. Rotor tip clearance flow 

4. Compressor bleed, stator leakage, and tur- 
bine cooling flow through inflow and out- 
flow leakage boundary conditions or source 
terms 


5. Stator button geometry 

6. Real gas effects with a linear relation 
between the ratio of specific heats and 
temperature 

7. Flow-path modifications directed by work 
of Hunter (ref. 4) to simulate cavities 

All turbomachinery components have been 
simulated at engine conditions. 

The HPC simulation in this full engine simula- 
tion builds on the HPC component simulation 
presented by Adamczyk (ref. 2). To match the 
thermodynamic cycle model of the HPC, 
APNASA was run by adjusting the inlet total pres- 
sure to match the exit corrected mass flow. The 
turbine simulation is identical to that presented by 
Turner etal. (ref. 10) except that the combustor 
profiles have been used as a boundary condition 
and the shaft speeds were set to the cycle values. 
The corrected and physical shaft rotational speeds 
were held constant at cycle values. 


Combustor Simulation 

The combustor has been simulated using the 
NCC explained by Ebrahimi et al. (ref. 3), Liu 
(ref. 6), and Ryder and McDivitt (ref. 8). This is a 
parallel-unstructured solver that uses a precondi- 
tioner to efficiently handle low Mach number 
flows. The Navier-Stokes flow solver is based on 
an explicit four-stage Runge-Kutta scheme. Turbu- 
lence closure is obtained via the standard k-E 
model with a high Reynolds number wall function 
or a nonlinear k-E model for swirling flows. The 
solver can be linked to any computer-aided design 
system via the Patran file system. CAD geometry 
can then be imported into the NCC system rela- 
tively easily. 

The NCC can be run with a gaseous fuel or by 
modeling the spray combustion process. Gaseous 
fuel with finite-rate chemistry has been used in the 
simulations thus far. 

Coupling 

The coupling between APNASA and NCC 
takes place at an interface plane. There are many 
issues to address when coupling codes that are 
very different, and even have different models for 
properties. Several key quantities must be 
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Fan, block 2 



Figure 2. — Five domains simulated in sequence inlet to exit. 


conserved from one code to the next. The mass 
flow must be conserved and is the most critical. 
The next most important quantity to conserve is 
the mass averaged total enthalpy. Since mass flow 
is conserved, this guarantees the total enthalpy 
flux is also conserved. For turbomachinery appli- 
cations, the next important conserved quantity is 
the angular momentum. After this, the total pres- 
sure is conserved, and finally the profiles of these 
quantities. Turner et al. (ref. 13) have presented 
details of this coupling. 

Figure 2 shows the five domains simulated. Each 
component had been mn near the cycle condition, and 
matched fairly well. The simulation is then mn se- 
quentially with the fan first, followed by the booster, 
HPC, combustor, and turbine. Scripts on the Ori- 
gin 3000 at NASA Ames have been written so the 
downstream domain boundary condition is set up and 
the downstream component is submitted to the batch 
queue automatically after the solution of the current 
domain has completed. 


Blockage Analysis 

Blockage is derived as a way to integrate 
overall three-dimensional effects in the tangential 
direction. It is often a quantity used as a catchall 
for analyzing experimental results when trying to 
match continuity and using measured flow angles, 
total temperature, and total and static pressure. 
Adamczyk (ref. 2) presented some spanwise inte- 
grated blockage calculations just downstream of 
blade row trailing edges for the GE90 compressor. 
Turner et al. (ref. 13) have presented the blockage 
for the GE90 turbine. The following analysis dif- 
fers from those definitions of blockage in that the 
mass-averaged total pressure has been used instead 
of the momentum-averaged total pressure. This 


was done because a mass average has somewhat 
more meaning in a turbine. Also the radial flow 
angle uses the area averaged axial and radial 
fluxes rather than finding the radial angle that 
would have zero normal momentum flux. This was 
done for convenience of the calculation. 

The following are defined at each spanwise 
grid-face of an axisymmetric grid by integrating 
the flow field: the mass-averaged values of total 

pressure ( P T ), total enthalpy ( h T ), and angular 

momentum (tVq ). The mass flow rate (m ) evalu- 
ated by integrating in the tangential direction and 
the area projections in the axial and radial direc- 
tion ( A z and A r ) are also needed. The following 
area averages are also used: radial momentum 
(pV r ), axial momentum (A z pV z ), and pressure 

( p ). The area averages are obtained using the 
weights and pointers needed for the closure as 
described by Kirtley (ref. 5). The mass-averaged 
quantities use area averages of fluxes. Quantities 
with an under-bar are derived. 


II 

£>1 

(1) 

<)) = tan 1 -BE 

Ip^J 

(2) 

v=^yI+a+vI 

( 3 ) 

V r = V z tan <]) 

( 4 ) 
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Figure 3. — Schematic of simulation. 


From equations (3) and (4), therefore 


y z =. 


y1-v£ 


1 + tan (j) 


The total enthalpy is defined such that 


h = h j 



(5) 


( 6 ) 


T and T t are calculated from h and h T using the 

linear function of y with temperature as described 
by Turner (ref. 10). An approximation to the isen- 
tropic relation using an average y yields 


P_ 

P T 



(y/(Y— i)) 


(7) 


Equations (5) to (7) and the relations of y on 

temperature are iterated because of the nonlinear 
enthalpy-temperature relation. The equation of 
state is 



The flow rate is 

™ — ^ aero^ metal P (YzAz ) (^) 


where b meta i is the metal blockage in the tangential 
direction. The term b aero is the aero blockage and 


has been introduced as an overall parameter to 
represent the three-dimensional effects. Solving 
for aero blockage yields: 


b 


aero 


ril 

b metal Pty z^ z + V r A r ) 


( 10 ) 


The derived absolute and relative flow angles can be 
calculated: 


3 

u. 

1 

II 

(ii) 

-i v e 
a = tan = 

(12) 

Ya 


a -i w e 

P = tan l = 

(13) 

V 


Results 



The simulation has fully coupled the fan, 
booster, HPC, combustor, and the turbines at the 
sea-level takeoff condition. Figure 3 is a schematic 
of the overall simulation result showing the axi- 
symmetric solution with an engine cross section 
superimposed. 

Figure 4 shows the geometry of the GE90 
combustor, where the dual annular fuel nozzle 
arrangement can be seen. The high-pressure tur- 
bine nozzles can also be seen in this figure, but are 
not part of this combustor simulation. The unstruc- 
tured grid used for the combustor simulation is 
shown in figure 5. Figure 6 shows contours of total 
pressure and total temperature at a plane cutting 
through the centerline of each fuel nozzle. The 
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Figure 4. — Geometry of the GE90 Combustor. 


Grid generated with FEMAP, adapted 
grid with NCC solution 



1 ,1 00,000 Tetrahedron 
24° periodic sector 


Figure 5. — Unstructured grid of the GE90 combustor 
used by NCC. 




Figure 6. — Combustor simulation. 
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Figure 7. — Full engine simulation comparison to cycle. Percent difference in total pressure (P), total 
temperature (T), and flow rate (W). 


TABLE I.— COMPUTER SIMULATION RUN TIMES. 


Component 

Number of 
blade rows 

Number of 
iterations 

Number of 
processors 

Wall-clock 

time 

Case from scratch (extrapolated) 

Number of 
iterations 

Number of 
processors 

Wall-clock 

time 

Fan 

3 

50 

16 

0:02:04 

8 000 

64 

1:22:40 

Booster 

8 

6 000 

256 

1:15:01 

6 000 

256 

1:15:01 

HPC 

21 

10 000 

512 

1:53:17 

10 000 

512 

1:53:17 

Combustor 


1 000 

256 

0:11:07 

31 000 

256 

3:53:00 

Turbine 

18 

10 000 

512 

1:56:36 

10 000 

512 

1:56:36 

Total 


5:18:05 


10:20:34 


dilution jets and their impact are clearly seen. Also 
seen is the substantial drop in total pressure across 
the combustor nozzles and liner. 

Based on extensive efforts to reduce runtime 
by exploiting parallel computing, the simulation 
has been run in just over 5 hr as shown in table I. 
The fan and combustor did not need to be run very 
long since they were already at the converged 
point. A solution from scratch with components 
near the cycle conditions would take 10 hr and 
20 min. 

A comparison of the simulation with the cycle 
is shown in figure 7. Percent difference in total 
pressure, total temperature, and flow rate are 
shown. The profile out of the booster has been 
adjusted to obtain the cycle pressure and tempera- 
ture at the HPC inlet due to the sensitivity of the 
full 21 -blade HPC simulation. The flow in the 


HPC is low as is the pressure ratio. Other discrep- 
ancies can also be seen for some of the other com- 
ponents. The problem is that the simulated 
components cannot duplicate the cycle conditions 
with the prescribed geometry (including the vari- 
able geometry settings). Work is underway to map 
the components as minimaps so that the cycle 
can be rerun and be self-consistent with the 
simulation. 

Figure 8 shows the power balance that has 
been calculated (it is normalized by the fan and 
booster work). The power from the fan is based on 
the enthalpy rise, whereas the other components 
are based on torque from pressure and skin fric- 
tion. The power to pump the cooling flows has not 
been subtracted from the HPT work. The work by 
the HPC has been under predicted which is consis- 
tent with the pressure rise. 
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Power from fan is based on enthalpy rise. 
Power from HPC, HPT and LPT are based 
on pressure and skin friction torque. 



Figure 8. — Power balance (power normalized by fan and booster). Power to pump cooling flows not subtracted 
from HPT. Work by HPC under-predicted.. 



Blockage 

1 1.000 

0.975 

0.950 

0.925 

0.900 

0.875 

0.850 


Figure 9. — Contours of blockage in the fan domain. 



Figure 10. — Contours of blockage in the booster domain. 


To the authors’ knowledge, this is the first 
complete simulation of an entire dual-spool engine 
at a given point. As mentioned, there is an on- 
going effort to make this simulation more consis- 
tent. These simulations present the state of the art 
for predicting the full engine aerodynamics. It can 
also be used for better component simulations, and 
to understand the flow fields better at realistic 
conditions. Blockage and the flow details it repre- 
sents will be presented as one method to look at 
the flow fields in the turbomachinery. 

Figure 9 shows contours of blockage for the 
fan domain that includes the fan, OGV and booster 
stator 1 . Clearly seen is the blockage due to the 
shock in the fan, the tip clearance vortex, and the 
secondary flows at the hub for both the fan and 
OGV. Contours of blockage in the booster are 
shown in figures 10 and 11. A separation near the 
casing of booster stator 2 is obvious from 
figure 1 1 . A radial profile of blockage and tangen- 
tial flow angles just downstream of booster sta- 
tor 2 are shown in figure 12. The tangential angle 
labeled as “derived” is based on equation (12), and 
the one shown as “based on momentum area aver- 
age” is based on the arc tangent of the area 



Figure 1 1 . — Enlarged view of Figure 1 0 showing blockage 
in the front part of the booster domain. 


average of the tangential and axial momentum. 
The difference in this angle is quite large where 
the blockage is large, and represents the impor- 
tance of the deterministic source terms in the aver- 
age passage model. Streamlines in the separated 
region of booster stator 2 are shown in figure 13. It 
is a comer separation that does not impact the loss 
immensely, but is important to the blockage and 
most likely to the stability of the machine. 
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Tangential flow angle, degrees 

Figure 1 2. — Radial profiles of blockage and tangential 
flow angle just downstream of the booster stator 2 
trailing edge. The derived angle is part of the blockage 
calculation. 



Figure 1 3. — Streamlines at casing of booster 
stator 2 showing corner separation. 



Blockage 



Figure 14. — Blockage contours for HPC, and zoomed 
region near stator 2. 


The blockage contours for the HPC are shown 
in figure 14, with a zoomed region near stator 2. 
Observe how the level starts at 1 at the inlet and 
goes to mean-line value of 0.88 by stage 10. This 
axial blockage distribution effects the matching in 
a highly loaded compressor such as the GE90. 
Figure 15 shows the radial profiles of blockage 
and tangential angles similar to figure 12, but 
behind the trailing edge of HPC stator 2. The aero 
blockage has a minimum value of 0.895 at 
70 percent span due to a separation on the suction 
surface of the stator from 50- to 90-percent span. 
The angles at the trailing edge are plotted in 


figure 16 that correspond to the radial profile of 
angles in figure 15. The contours show under- 
turning at the hub and casing. The separation is not 
very conspicuous in this plot due to the large ra- 
dial exit angle variation. However, it is useful to 
see the extent of the tangential variation of angle. 

Figure 17 shows the blockage in the HPC 
domain in the vicinity of rotor 6. Fow values of 
blockage can be seen near the hub and tip of 
rotor 6. Figure 18 shows the cause of this block- 
age: a separation on the suction surface near the 
hub and the tip clearance vortex. Figure 19 shows 
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Figure 1 5. — Radial profiles of blockage and tangential 
flow angle just downstream of the HPC stator 2 
trailing edge. 


the radial profiles of blockage and absolute tan- 
gential flow angle. Even at midspan, where the 
flow is clean, the blockage has a value of 0.94. 
The cumulative multistage effects, through the 
deterministic sources have contributed to this low 
value. An angle difference of 2° to 3° can be seen 
between the blockage derived quantity and that 
based on the area averaged momentum. 

Figures 20 and 21 show the contours of block- 
age below unity for the HPT and LPT, respec- 
tively. There are values above unity that exist from 
just upstream of the throat to downstream of the 
uncovered region for most of the blade rows, es- 
pecially the nozzles. It is not fully understood why 
the blockage would be greater than one, and 



Figure 1 6. — Absolute flow angle at the 
trailing edge of stator 2 of the HPC. 

Angles vary from 1 5 (blue) to 30 (red) 
degrees. 

Blockage 
- 1 .000 

■ 0.975 

0.950 

0.925 

0.900 

0.875 

* 0.850 

Figure 17. — Blockage contours in the HPC domain near 
rotor 6. 

is currently being explored. It is consistent with 
turbine flow coefficients being greater than unity 
under some conditions. These contours still show 
areas to explore. One of these is a hub secondary 
flow vortex in LPT rotor 5 shown in figure 21 and 
through vorticity contours and streamlines in 
figure 22. 

In highly loaded turbomachinery, blockage plays 
a critical role in the radial distribution of flow and 
being able to pass the flow, especially at high Mach 


/—HPC rotor 6 
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Figure 18. — Streamlines in HPC rotor 6 
showing tip vortex and hub corner 
separation. 


numbers. As demonstrated, aerodynamic blockage 
demonstrates the location of shocks, secondary 
flows, separated regions, large wakes, and multi- 
stage interaction effects. These multistage effects 
come from the deterministic sources that are part 
of the average passage approach. A mixing-plane 
approach would produce a loss, and a discontinu- 
ity at the interface, and not a blockage. The aero- 
dynamic effect is not the same. Often a 
compressor rotor and stator are analyzed as iso- 
lated blade rows. It is very difficult to get the ef- 
fects of blockage that are implicit in the average 
passage solution, except by running at a lower 
flow and lower speed and modifying the boundary 
conditions. However, this does not fully mimic the 
total effect of blockage since the machine runs 
differently at higher speeds and the variation of 
blockage is not uniform spanwise as shown for a 
few of the blade rows. Flow angles for the stator 
are also very difficult to get right when the level of 
blockage is not correct. Turbines are also affected 
by blockage such as turbine throat setting and 
turbine thrust balance. These can be predicted 
better when blockage is addressed in the design 
process. 
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Figure 19. — Radial profiles of blockage and tangential 
flow angle just downstream of the HPC rotor 6 trailing 
edge. 


Conclusions and Future Work 

A successful simulation of the three- 
dimensional flow in the primary flow path of the 
GE90-94B high bypass ratio turbofan engine has 
been achieved. It would take less than 1 1 hr of 
wall clock time if starting from scratch through the 
exploitation of parallel processing. The engine 
simulation matches the engine thermodynamic 
cycle at a sea-level takeoff condition. Compari- 
sons with the cycle are presented. Details of the 
turbomachinery solution using a blockage analysis 
have also been presented. Aerodynamic blockage 
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Figure 20. — Blockage contours in the HPT. 
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Figure 21 . — Blockage contours in the LPT. 


demonstrates the location of shocks, secondary 
flows, separated regions, large wakes, and multi- 
stage interaction effects. 

Future work will couple the simulation to the 
cycle analysis using minimaps of the simulated 
components. The cycle will be used to balance the 
simulation and a front-to-back simulation will be 
performed. If necessary, this process will be iter- 
ated. This should eliminate any flow discrepancies 
between components in pressure, temperature and 
flow rate seen in the current simulation. 
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